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Abstract 

We study the energy flow due to the motion of topological solitons in nonlinear extended systems 
in the presence of damping and driving. The total field momentum contribution to the energy flux, 
which reduces the soliton motion to that of a point particle, is insufficient. We identify an additional 
exchange energy flux channel mediated by the spatial and temporal inhomogeneity of the system 
state. In the well-known case of a DC external force the corresponding exchange current is shown 
to be small but non-zero. For the case of AC driving forces, which lead to a soliton ratchet, the 
exchange energy flux mediates the complete energy flow of the system. We also consider the case 
of combination of AC and DC external forces, as well as spatial discretization effects. 

PACS numbers: 05.45.Yv, 05.60.-k 
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I. INTRODUCTION 



Solitary waves or solitons represent one of the striking and famous aspects of nonlin- 
ear phenomena in spatially extended systems. They appear as specific types of localized 
solutions of various nonlinear partial differential equations and possess several important 
properties ]] • These properties include an exponential localization of energy remaining un- 
changed during soliton propagation, and an elastic scattering of solitons. The latter one, 
i.e. the fact that two solitary waves maintain their identity after collisions, inspired the in- 
troduction of the term soliton |2j, which emphasizes the analogy between these objects and 
particles. Later it has been realized, however, that an essential condition for the particle-like 
behavior of solitons is the integrability of the corresponding models. But the effect of dy- 
namical localization of energy was found to be much more general, and the notion of soliton 
has been extended nowadays to non-integrable models as well. An analogy between such 
"generalized" solitons and particles is less obvious. Not only the scattering process is gen- 
erally different from the elastic scattering of particles, but several other important physical 
effects may appear for solitary waves in non-integrable models. An important feature essen- 
tially breaking the particle approach to soliton dynamics is connected with the possibility 
of excitation of the so-called shape 3i] or internal [4] modes: in contrast to a rigid body, a 
soliton may perform shape oscillations while propagating. The mechanism of possible energy 
exchange between the soliton kinetic energy, associated with its translational motion, and 
the internal energy of shape modes, can drastically change the soliton dynamics and leads 
in some cases to resonances during soliton-soliton interactions jjj. 

Another important aspect of soliton dynamics in non-integrable systems, which we ad- 
dress in this paper, concerns the problem of energy transport associated with a directed 
soliton propagation. Since a soliton carries some (finite) non-zero energy, it is naturally to 
assume that its translational motion would result in a non-zero energy current in the system. 
One can introduce a total field momentum, in analogy to a particle mechanical momentum, 
which precisely corresponds to the energy current when considering solitons in integrable 
systems [if. However, an attempt to extend this momentum-based approach to dissipative 
solitons may lead to confusing and misleading results. As an example, we mention here 
the so-called soliton ratchets: the effect of a unidirectional motion of a topological soliton 
(kink) under the influence of a spatially homogeneous and time-periodic external bias with 
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zero mean 



B B fl 3 S 10, II, 12, Q- As an extension of single particle ratchets Q, 
at a first glance soliton ratchets demonstrate at least qualitatively a nice analogy between 
solitons and particles. Indeed, applying the same symmetry analysis as for particle ratchets 
, one observes a directed motion of kinks not only in numerical simulations, but also in 
experiments with annular Josephson junctions However, a surprising and puzzling 

result appears if one estimates the energy current for the soliton ratchet case: simple ana- 
lytical calculations tell that the averaged value of the total field momentum is strictly equal 



to zero 
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17| . An analogous effect of non-topological soliton motion with zero total field 



momentum has been reported also for a nonlinear Schrodinger equation with parametric 
driving in the presence of dissipation [18]. These results show, that the total field momen- 
tum is no longer an adequate quantity to describe the energy flow due to soliton dynamics 
in non-integrable damped and driven systems. 

In this paper we study the energy flow due to soliton dynamics in the framework of the sin- 
Gordon model with a spatially homogeneous external driving and damping. We show, that 
in general the energy current associated with a translational motion of the topological soliton 
(kink) consists of two contributions. The abovementioned total field momentum constitutes 
only one component of the total current - the internal current. Another path for energy 
transport in the system is mediated by an inhomogeneous in time and space energy exchange 
between the soliton and the external degrees of freedom (driving force and damping). The 
corresponding exchange current has no analogy within the particle description of a moving 
soliton. Its existence is solely due to the spatial extent of the soliton. Even for the case of a 
time-independent driving, when shape modes are not excited, the exchange current is small 
but nonzero. It is drastically enhanced by the excitation of shape modes via a time-periodic 
driving or a spatial discretization of the system. 

The paper is organized as follows. In Sec. |H] we introduce the driven and damped sin- 
Gordon model and the basic definitions of the energy balance equation and the total field 
momentum (internal energy current). In Sec. IIHI the case of soliton motion under the 
influence of a DC external force is discussed. We demonstrate, that the internal energy 
current does not account for the total energy flow in the system. We introduce the notion 
of the exchange current, which completes the energy flow balance. A generalization of our 
approach to the case of soliton motion induced by AC driving forces (soliton ratchet) is given 
in Sec. EXl The effects of discretization are discussed in SecEl Finally, Sec. EH concludes 
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the paper. 



II. THE MODEL 



We consider a topological soliton (kink) motion in the driven and damped sin-Gordon 
model |l9j, widely used in the field of soliton ratchets , which describes 

;he dynamics of the superconducting phase difference across the annular Josephson junction 
3: 



<Pu ~ (fixx = -ot(pt - sin ip + E(t) . 



(1) 



This equation describes the evolution of a scalar field ip in space and time. Here and in 
what follows subscripts x and t denote partial derivatives with respect to the corresponding 
variables. The parameter a regulates the strength of damping in the system, and E(t) is an 
external driving force. In addition, we impose the kink-bearing periodic boundary condition: 



ip(x + L,t) — ip(x, t) + Q , (p t (x + L, t) — (p t (x, t) 

where Q = 2tt is the topological charge and L is the system size. 
The field energy density is given by 

p[<p(x, t)} = p{x, t) = + <pl) + 1 - cos(^) . 
Its dynamics is governed by the following energy balance equation: 



(2) 



(3) 



Pt 



-jl - (*Pt + E{t)ip t . 



(4) 



Here we introduce the internal energy current J 7 (t) and its density j ! (x,t): 



At) 



fdx , j 7 (x,t) = -p 3 



Pt 
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(5) 

By differenti- 



J 1 is also known in the literature as the total field momentum 
ating Eq. (EJ) and using Eq. (JTJ) together with the boundary condition ((2j) it is straightforward 
to show |l 7L I20I 21] that the internal current J 1 satisfies the following ordinary differential 
equation: 



j/(t) = -aJ\t) - QE(t) . 



(6) 
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In the absence of an external force E(t) = and dissipation a = (i.e. in the non- 
perturbed integrable sin-Gordon model), Eq. (£[]) with the above boundary conditions (J2J) 
supports the well-known kink solution which in the limit of an infinite system size 
L — > oo takes the form 

' x-Vt 



(f kmk (x, t) = 4 arctan ^ exp 



(7) 



Vi -v 2 

where V is the kink velocity, \V\ < 1. According to the energy balance equation (@J), in 
this case the internal energy current J 1 is the only possible pathway to mediate an energy 
transport in the system. Using the expression (JJJ) for the kink solution and the definitions 
of the internal energy current J 1 (J3J) and energy density p (jSJ), it follows, that the energy 
current associated with the kink motion can be obtained in full analogy with a moving point 
particle: 

J 1 = VW kink . (8) 
Here the kink energy W kmk is obtained from the energy density p(<p): 

W kmk = £ p[ V kmk ]dx , W k Zt = J^ys • (9) 

In the general case, when both E and a are non-zero, a moving kink solution of Eq. (JT|) 
can persist as an attractor of the system with all the parameters (including the kink velocity) 
determined by the choice of E, a and L. Our goal is to compute the energy current generated 
by such a moving kink. As will be shown below, an inhomogeneous exchange of energy 
between the kink and the external degrees of freedom (force and dissipation), opens an 
additional path for energy transport in the system, which, together with the above internal 
current, constitute the full energy current balance generated by the moving kink. 



III. CONSTANT DRIVING FORCE 

Let us start with the seemingly simple case of the perturbed sin-Gordon equation (JTJ) - 
the case of a constant driving force E(t) — E = Const. A kink moves in this case with a 
constant velocity V (defined by a choice of the force strength E, as well as the damping 
constant a and the system length L) j _so that the corresponding attractor solution depends 
only on a single variable £ = x — Vt : 

(p(x,t)=^(x-Vt)=^)' ( 10 ) 
5 



Due to the time homogeneity the total energy of the system W, 

W = f L p[tp(x,t)]dx , (11) 





remains constant in time on the attractor solution. For a point particle the current balance 
would not change, since the energy loss due to the damping would be compensated by the 
external dc field. For the spatially extended kink this will not be the case anymore. 

Let us assume, that the system size L is chosen to be large enough, so that far away from 
the kink center the spatial field distribution asymptotically approaches the homogeneous 
ground state: — >• 0) — > <p v , —>■ L) — ► <p v + Q. Note, that the external force E 
shifts the degenerate ground state of the system from ipQ = 2irm,m = 0,±1,±2,... to 
= ^ + <p v (E), <p v (E) = arcsin(E). 

We denote by w(x, t) the amount of energy stored on the stripe [0,x] at time t. We 
further assume that x = is far from the kink center so that the energy current density 
vanishes there. Then we can define the local energy current density j(x,t) as 

j(x,t) = = -£ Pt ( x >,t)dx' . (12) 

Using the attractor property (|1U|). we obtain: 

j(x-Vt)=j(£) = V\p(t)-p(p)] , (13) 



and the total energy current J is thus given by 



Jo Jo 



(14) 



The quantity in the brackets in the r.h.s. of Eq. (J 14)) is the difference between the system 
energy with and without the kink. Indeed, in the presence of external force and damping, the 
ground state energy has shifted from zero to Wq = L ■ p((p v ) = L ■ p{x = 0). This difference 
is precisely the kink energy: W k = W — Wo- Thus, we arrive to the simple relation between 
the total energy current J, kink velocity V and kink energy W k : 

J = VW k . (15) 

As expected, this expression corresponds to that of a point particle moving with velocity V 
and carrying energy W k . However, opposite to the integrable case E = a = 0, the internal 
current J 1 does not coincide with the total current J: J 1 ^ J. Indeed, using the definition 



of the total current J given in Eqs. (jl2|) . (jl4j) and the energy balance equation (0J), we arrive 
at 



J = J 1 + J E , (16) 



where the exchange current J E [l6j] has been introduced 

"L px 



J E = V [ L dx f di \aV (V/(0) 2 + ■ (17) 

Jo Jo J 

Following the definition of the internal current J 1 in Eq. (jSJ), we obtain 

J 1 = V [ L ' {ip'{x)f dx , (18) 



so that all components of the current balance equation (fT^j) are expressed through the field 
ip. Below we will calculate numerically the moving kink attractor solution ip(£) in order 
to estimate all the three currents for different values of the model parameters. Note, that 
according to Eq. (jUJ) one can obtain the exact value for the internal current J 1 on the 
attractor solution in an independent way: 

J exact = • (19) 

a 

We will use this expression in order to control the accuracy of our numerical schemes. 

The fact that the moving kink solution <ft(x, t) depends only on the single variable £ = 
x — Vt allows to reduce the original field equation (0) to the well studied driven and damped 
pendulum problem. Indeed, substituting ansatz (fTU|) into Eq.(^), we arrive to the following 
ODE for the function 

(1 - V 2 )ip" = -aV$ + sin(V>) - E, (20) 

which describes the evolution of an effective pendulum with the momentum of inertia M = 
(1 — V 2 ). The damping constant 7 = aV depends on the kink velocity V. The external 
force E plays the role of an applied constant torque to the pendulum. By introducing the 
effective time 



Vi- v 2 ' 



(21) 



Eq . (|2(J|) can be finally re- written for the function y(r) = ip(^/\/l — V 2 ): 

y = -Gy + sm(y) - E, (22) 



with the damping constant G defined through the original kink parameters: 

G{E, L) — (23, 
1 — V 2 (E, L) 



Here we emphasize the implicit dependence of G on the strength of the external force E and 
the system size L, since the kink velocity V is unambiguously defined by these parameters. 

Thus, a moving kink solution of the original problem Q with the boundary con- 
ditions (j2J) corresponds to a limit cycle solution y(r) of the driven and damped pendulum 
()22j) . On this limit cycle the pendulum performs periodic rotations with the period T defined 
through the original system length L and the kink velocity V as: 

T = - L (24) 

Note also, that the limit of an infinite system size L corresponds to a limit cycle solution of 
the pendulum with an infinitely large period T. 

The existence and the characteristics of the limit cycle in Ea. (l2*2l depend on the choice 

nnn , , 

of values of G and E 24, 25, 26J. For each given value of E from the interval \E\ < 1 
there is a critical value of G = G cr (E), below which the limit cycle co-exists together with 
a stable static pendulum state (the latter corresponds to the homogeneous ground state in 
the kink tails). The period of the limit cycle tends to infinity as G approaches G cr from 
below, while above the critical value of G the limit cycle disappears. As a consequence, the 
value G cr (E) corresponds to the limit of an infinite size L in the original problem (0) with 
the corresponding drive E, while any lower value G < G cr is associated with a moving kink 
in a finite-size system. For small values of \E\ the critical value of the effective damping can 
be approximated as 2(| 

7lE 

G cr « — , \E\ < 1. (25) 

By using relation ()23|) one can derive the corresponding approximation for the kink velocity 
in an infinite-size system at low power of the external drive: 

y» - — L= LE|<1. (26) 

To find the limit cycle solution of ()22|) and the corresponding value of G cr at arbitrary 
values of E we use the Newton method in phase space (y,y), rewriting Eq.(|22} as j3]: 

z(y) = y, (27) 

^ = -G + ^^, (28) 
ay z 
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FIG. 1: Kink velocity as a function of the dc force E for a = 0.2. Open circles: numerical solution 
of (|22|) . see text for details, solid line: approximative solution (|26|). Inset: critical value G cr from 
(|23|) for a = 0.2. Open circles: numerical solution (see text for details), solid line: approximative 
solution <|2"5|). 



and computing the map 



{^(2/o), 2/o} 
2/o 



{z(y + 2n),y + 2n} 
arcsin(i?) . 



(29) 
(30) 



For any z = z(yo) a unique periodic solution of ()27I28|) satisfying condition z(yo + 2n) = z 
can be found by tuning the parameter G, provided that the sign of Zq is chosen in accordance 
with the sign of E. The critical value G cr is obtained by taking zq = 0. Note, that the value 
of zq is proportional to the value of = 0,L). For any finite system size L it is nonzero 
due to exponentially decaying kink tails. 

In Fig. ^ the computed values of G cr , as well as the corresponding values of the kink 
velocity V, are shown with circles. In fact, the approximation for the kink velocity ()26|) 
works reasonably well even at values of E close to unity, while the critical value G cr deviates 
stronger from the linear approximation (|25|) at values of the driving force close to E — 1. 

Once the limit cycle solution y{r) of (|2"2"j) is found, the values of the total (fPfj) . exchange 
(fTTjl and internal (|T%|) currents for the corresponding moving kink solution can be obtained: 

aG 



J 1 
J 



a 2 + G 2 
2G „ 



(D-F) 



a 



Ga 2 + 2G 2 
a a 2 + G 2 



F + 



aG 



a* + G 2 



D 



(31) 
(32) 
(33) 



= J E + J 1 , 
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where the functions F and D, are related to the mean pendulum kinetic and potential 
energies: 

F(E,L) = [ dr V -, (34) 
o 2 

T 







D(E, L) = / dr Vl-E 2 - cos(y) . (35) 



Here the implicit dependence on the original system length L comes through the pendulum 
period T. 

Therefore, the three functions D (JH3J), F (jH3| and G uniquely defined for a limit cycle 
of (|22|) for a given value of the external field E and 2/(0), determine all the characteristics 
of the corresponding moving kink problem: its spatial profile, velocity, energy, internal, 
exchange and total currents. In Table H] we list the values of all three currents calculated 
for the damping constant a = 0.2 and different strength of the external driving force E 
by means of the pendulum approach, corresponding to the limit of L — > oo in the original 
system (JIJ. They perfectly agree with the results of direct numerical integration of Eq. {1}. 
In order to control the accuracy of our numerical scheme, we estimate the error of the 
calculated internal current |Z\ J" 7 1 , for which we know the exact result Ji xact ©• We also 
indicate in Table |J the accuracy A J at which the current ballance equation (fTT)|) is fulfilled: 
A J = J — J 1 — J E , and the relative strength of the exchange current \ J E / J\. 

The surprising result of the presented analysis is that for the case of a constant driving 
force, when the kink is known to behave similar to a point particle without any internal 
degrees of freedom being excited, the exchange current is nonzero! It is relatively weak, but 
both its absolute value and its relative strength are significantly larger than the estimated 
numerical error and increase with increasing driving force E. 

Since for a constant driving force all currents are time- independent, we can easily define 
the densities of internal and exchange currents 

ji{x) = V^\xf , (36) 

f(x) = V f di [aV (V/(£)) 2 + E^'(0] , (37) 
Jo 1 J 

using Eqs. (fT£J) and (fTTj) . The corresponding profiles are plotted in Fig. El for the case 

a = 0.2, E — 0.2. First it is evident from (j37jl . that the exchange current density is nonzero, 

similar to the internal current density (|36|). This is due to the finite spatial extent of the 

kink, which causes a corresponding spatially (and thus temporally as well) inhomogeneous 
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TABLE I: The values of the internal, exchange and total currents for the damping constant a = 0.2 
and an infinite system size L — > oo, calculated within the pendulum approach . A J 1 and A J are 
the errors of the calculated internal current and the accuracy at which the current balance relation 
(|16|) is fulfilled, respectively (see the main body text for details). 

energy exchange between the field if, the external forcing E and the damping term. A total 
vanishing of J E could only hold due to some spatial symmetry of j (x) which in turn has to 
be generated by some symmetry of the kink profile. However the kink profile in the presence 
of driving force and damping is asymmetric. That is the main reason for the nonvanishing 
exchange current contribution. 

Finally, we note, that the pendulum approach is valid only for the case of \E\ < 1. Indeed, 
in the case of \E\ > 1 the limit cycle of (j22j) persist, but the stable fixed point disappears. 
This means that one no longer has a homogeneous ground state in the kink tails. In this case 
the ansatz (jl()j) is not valid and the pendulum approach becomes inappropriate. Similarly, 
when \E\ > 1 the proposed scheme for the calculation of the total and exchange currents, 
Eqs. ()14j) and ([17)1 . breaks, since the homogeneous ground state in the kink tails is no longer 
supported by the starting equations. 
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FIG. 2: The densities of the scaled internal current, 0.02 j 1 (x) (thin line), and of the exchange 
current, j E {x) (bold line), for the E = 0.2, a = 0.2 and L = 500. 

IV. TIME-PERIODIC DRIVING FORCE: SOLITON RATCHETS 



In this section we proceed to the case of a time-periodic driving force E(t + T) = E(t) 
with zero mean value. For the case of an AC driving force, which possesses the time-shift 
symmetry 



E(t) = -E(t + T/2) 



(38) 



the combined symmetry transformation 

T 

x — > —x, ip — > —<f + Q, t — >• t + — 
leaves Eq. (0) invariant while changing the sign of the kink velocity V denned e.g. as 



1 r L 

V(t) = — / xiptxdx . 
Q Jo 



(39) 



(40) 



If only one kink attractor persists, then on this attractor the kink will not move, but oscillate 
at the best (IJ Q] . This is the typical situation observed in various numerical simulations. In 
the unlikely case that a kink attractor persists on which the kink is moving with a nonzero 
average velocity, then by symmetry another attractor exists as well, on which the kink 
moves with the opposite velocity. However, the above time-shift symmetry (J38|) is generally 
removed by choosing the bi-harmonic driving force 



E(t) = E 1 cos[cu(t - t )} + E 2 cos[2cu(t - t ) + ©] , w 



2tt 



< 6 < 2tt . 



(41) 



12 




(a) 



(b) 



260 



255 



X 250 



245 



240 



260 



255 I 



X 250 



245 



240 




FIG. 3: Space-time evolution of the soliton ratchet with the E\ = E\ = 0.2, uj = 0.1, = and 
a = 0.2, L = 500: (a) Contour plot of the energy density ©; (b) Contour plot of the function 
4>{x,t) defined in Eq. (j5U|) . 



As a consequence, the kink on the attractor solution will in general have a nonzero velocity 
(provided no other hidden symmetries of the system persist). This results in a unidirectional 



(in average) motion of t 
given in Eq. (j^T| 
numerical simulations 




le kink under the influence of zero-mean periodical driving force 



1^ . The corresponding kink ratchet effect was observed both in 



and in experiments with an annular Josephson junction 12]. In 
Fig. ED^a) the characteristic evolution of the energy density (JHJ) in the kink ratchet dynamics 
with the bi-harmonic driving force (j41j) is shown. A gradual drift of the energy density peak, 
associated with the kink center, is clearly observed. 
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The unidirectional motion of the kink should result in the appearance of a non-zero 
averaged energy current. Simultaneously, as follows from Eq.ljfi]). for any zero-mean E(t), 



the time average value of the internal current J 1 should be zero on the attractor 17ll20ll21|. 
Therefore, the internal current does not contribute to the energy transport associated with 
the kink motion. This is impossible from the point of view of a point particle, since then 
the internal current J 1 corresponds to an effective particle momentum jigL 2(1 21, 22 1, and 
the kink will perform a momentumless motion in the ratchet case. This circumstance raised 
discussions of whether the observed rectification of the kink motion is due to additional 



white noise or discretization effects (see e.g. 13] )• However, such a momentumless motion 
does not contradict the energy balance principles, since there exists an alternative path for 
the energy flow in the system - the exchange current J E , as we demonstrated in the previous 
section. In the case of bi-harmonic driving force all the energy is transported through the 
exchange current, so that the time- averaged current balance equation (|16|) reads: 

J = J E , J 1 = . (42) 

Here we use calligraphic letters to indicate time-averaged quantities: J = (1/T) 2 • 
Jo dt[£ T dtxJfa) (details of the averaging procedure will be discussed below, see also 
Ref. Ilj). 

Our goal now is to estimate the energy current J generated by a moving kink. Taking the 
system size L much larger than the characteristic oscillation distance L p which the moving 
kink covers during a single period of the external force, we may assume that far from the 
kink center the field distribution in space is homogeneous. We separate the field variable 
(p(x,t) into a localized kink part, ip k (x,t), where ip k (x — *> 0,t) — 0; ip k (x — > L, t) = Q, and 
a background (vacuum) part (p v (t) which depends only on time |2l|: 

<p(x,t) = V k (x,t) + V v (t) . (43) 

The vacuum part must satisfy Eq.Q in the absence of the kink, i.e. when Q = 0. Therefore 
it does not contribute to any energy transport [sj]. 

On the attractor the dynamics of the system (|TJ is given by 

V k (x, t + T) = if k (x - VT, t) , <f v {t + T) = <f v (t) , (44) 

where V is the averaged kink velocity. Note that all integral system characteristics, such 
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as the total energy of the system, the kink velocity, the energy current, etc. are periodic 
functions of time with the period T. 

Since the total system energy W is not conserved, we can not introduce instantaneous 
energy currents. However, since W is a time-periodic function on the kink attractor solution, 
we estimate currents by considering energy density distribution changes over the period of 
attractor, when the total system energy is restored. Assuming that x = corresponds to a 
kink tail point, the amount of energy Aw(x,t): 

Aw(x, t) = f X \p(x', t) - p(x', t + T)]dx' = 
Jo 

- [ X dx' [ p t (x',t')dt' (45) 
Jo Jt 

which the part of the system [0, x] looses or gains during one period T has to be transported 
through the point x by the energy current j(x,t) = Aw(x,t)/T. Both quantities, j(x,t) 
and Aw(x,t), are obviously periodically dependent on the reference time t. However, it 
is missleading to associate them with any instantaneous characteristics of the system, as 
Aw(x,t) is defined on the timescale of the attractor period T. In order to eliminate the 
dependence on reference time t we average both quantities over the period T. Finally, the 
total averaged energy current in the system is obtained by the additional integration over 
the space: 

J = )-\ L dx [ 7 dtj{x,t) , j{x,t) = ^Aw(x,t) . (46) 
1 Jo Jo 1 

Using Eqs. (|44j) and the multiplicative integration rule 23| we finally obtain (see 
Appendix for details): 

J = V( [ L p{v(x,t)}dx - Lp[^(t)}) T (47) 
Jo 

with (...)t = 1/T Jq ...dt. As before, the difference between the total system energy W and 
the background energy Wq, averaged in the r.h.s. of Eq. ()47|) . denotes the kink energy W k . 
Therefore, the general result for the total energy current, given in Eq. (jl4J) . still holds in 
this case, if using averaged values. 

Using the results ()45|46j) we can now derive an expression for the averaged exchange 
current J E from the energy balance equation (0J): 

J E = - j dt [ dxf{x,t) , (48) 
1 Jo Jo 
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FIG. 4: Numerically computed average values of the exchange current J E (lines) and total current 
J (circles) as a function of the driving force parameter for a = 0.2 (solid line, filled circles) and 
a = 0.05 (dashed line, open circles). Other parameters are the same as in Fig. El 



Note, that similar to previously defined j(x,t) in Eq. ()47|). the quantity j E (x,t) can not 
be treated as an instantaneous current density value. 

Again, using the multiplicative integration rule j23], Eqs. (|4~K|) - (|oT)j) can be reduced to the 
more compact form (see Appendix for details): 



where (jf{t) = <f>[(f v (t)]. 

The appearance of a non-zero exchange current in the system is mediated by a spatially 
and temporally inhomogeneous energy exchange between the moving kink and the external 
degrees of freedom. The information about the energy exchange process is contained in the 
function <j>(x, t) (foTIj) introduced above. In Fig|3fb) the space-time evolution of <fi is plotted. 
The energy is exchanged and transported in a cyclic way: first the kink absorbs energy in 
its rare tail, then it releases energy in its front, then it absorbs energy in its front and finally 
releases energy in the rare tail. 

The numerical simulations of the kink dynamics (P) confirm that the balance between 
averaged values of the total and exchange currents holds, see Fig. |H The absolute values of 
the exchange current increased, as compared to the case of a constant driving force (cf. with 




(49) 
(50) 




(51) 
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Table It may be due to the excitation of kink shape modes, which are clearly observed in 
Fig. the kink is much more compressed when moving opposite to its average propagation 
direction as compared to the times when it moves in the same direction. 

Note that the rectification effect vanishes at a certain value of the phase mismatch 
between the two components of the driving force. This value depends on the damping 
constant a. That is due to another symmetry property of the system in the Hamiltonian 
limit a — > 0. Indeed, in that limit the system possesses the time-reversal symmetry t — ► —t. 
That symmetry operation changes the sign of the kink velocity ()40|). provided the driving 
force is taken to be symmetric in time: E(—t) = E(t). Therefore, in the Hamiltonian limit 
the rectification effect should disappear at 6 = 0, %. In the underdamped regime «Cl this 
effect persists at some value of 6 close to 6 = 0, 7r 

If a DC component is added to the bi-harmonic driving force ([41)1 . a careful tuning of 
the parameters can lead to an exact cancellation of both DC and AC force components, so 
that V = 0. Then the averaged total current J (14 7|) is exactly zero. However, according 
to Eq. (fTTJj) . the internal current is J 1 = —QE stop /a, where E stop is the DC component 
of the driving force. This implies an exact balance between the two nonvanishing current 
components: 

J = 0, J E = -J I . (52) 

For example, taking parameter values from Fig. EJ we found E stop pa 0.01705 and J E = 
—J 1 pa 0.536. Note also, that for any < E < E stop the sign of the total momentum of the 
system J 1 is opposite to the sign of the kink velocity V, and the internal current is pumping 
energy in a direction opposite to the kink motion. 



V. EFFECTS OF SPATIAL DISCRETIZATION 



Up to now we discussed the energy flow due to kink motion in a spatially extended con- 
tinuous system. It is known, that the spatial discreteness of a system can drastically change 
the soliton dynamics j^J. By breaking translational invariance of the system, discreteness 
induces a Peierls-Nabarro potential which strongly influences the translational motion of 



the kink 



. It is also known, that additional internal modes of the kink may appear in 



discrete systems |28j]. Their existence can essentially modify energy exchange mechanisms 
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between the moving kink and external driving forces and damping. As a consequence, the 
relative contributions of the internal and exchange currents to the total energy current can 
be strongly changed in the discrete case both for DC and AC driving forces. 

Let us consider the discrete version of Eq. (JTJ also known in literature as the damped 
and driven Frenkel-Kontorova chain 2^|: 

(p n - C 2 (<f n +l + (Pn-1 - tyn) = -OKpn ~ Sin(p n + E(t) . (53) 

It has been used in solid state theory for modeling dislocation dynamics in crystals. This 
equation also describes fluxon dynamics in one-dimensional arrays of coupled Josephson 
junctions j^]. It is important to note, that Eq. ()53|) is often used for numerical simulations 
of the original continuous sin-Gordon system (fl]). where C ~ 1/h is the inverse mesh size 
or the lattice spacing (i.e. C — > oo corresponds to the continuous limit). In this respect it 
is also important to investigate possible effects of discretization which might appear during 
the numerical calculation of energy currents of an underlying spatially continuous system. 

In order to relate the above results for the continuous sin-Gordon model (0) to the 
discrete system flSHj) . we use the transformation x — > n/C, <f(x,t) — > (p(n/C,t) = <p n (t). 
Thus, all spatial derivatives should be replaced by finite differences weighted by the factor 
1/h = C, while all spatial integrals should be replaced by discrete sums over lattice sites 
with the weight factor 1/C. With that we obtain the following discretized versions for the 
internal (0) and total ()47|) energy currents: 

j'^^Eln, (54) 

J = V d J2[Pn-Np v ], (55) 

n 

where Vd is the kink velocity calculated in units of lattice site differences per time. The local 
internal current and energy density p n are yet to be defined. 

We note, that the Frenkel-Kontorova model corresponds to the Hamiltonian 

H = Pn = + ~4~ [( < ^«+ 1 ~ V 9 ".) 2 + (<Pn ~ <Pn-i) 2 \ + 1 - cos(y? n ) . j (56) 

Following essentially the same steps as in Sec. |H1 we can derive the discrete version of 
the energy balance equation (JU) for the discrete energy density p n (energy of an individual 
lattice site): 



C (in" Jn-l) + (57) 



dpn 

dt 
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FIG. 5: The internal (filled circles), the exchange (triangles) and the total (squares) energy currents 
as functions of the discretization parameter C 2 . (a) and (b): the case of the DC external force 
E = 0.2; (c) and (d): the case of the biharmonic driving force (|41[) with E\ = E% = 0.2 and 
G = 0. a = 0.2 for all cases. The dash-dotted line in (b) indicates the exchange current value, 
J E ~ —0.0320128, calculated for the continuum limit within the pendulum approach (see Table 
1.). The dashed line in (c) corresponds to the power-law asymptotic J ~ C~ 2 . 



c 

J'n = 'J + <Pn) (Vn+1 ~ Vn) ■ 



(58) 



The internal current density (also known as the local heat flux j3l|) is the discrete analogue 
of the internal current density (jHJ) of the continuum case. 

In contrast to the continuum limit, it is not straightforward to define a discrete version 
of the exchange current on the basis of the last two terms in the r.h.s. of the balance 
equation (foTj). However, the exchange current can be computed via the difference between 
the total and internal currents, J E = J — J 1 . 

In Fig. the dependencies of the computed energy currents on the discretization param- 
eter C are plotted for the cases of DC and AC external driving forces. Increasing C, in both 
cases all the currents clearly converge to the corresponding values obtained in the continuum 
limit: in the DC case the dominant energy pathway is realized through the internal current, 
while in the AC case the only remaining path for energy transport is mediated by the ex- 
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change current. In the latter case the internal current scales down to zero as J 1 ~ C 2 for 
C2 > 10, w hlch foUows fr0I „ the app—on of integrals by d.serete 8Um8 Q. Spatia, 
discretization can drastically change the ratio between the two currents when C 2 < 10 both 



in the case of a constant external bias, and in the case of a soliton ratchet 
these values of C the Peierls-Nabarro potential is still exponentially small 



3 



Since for 



, the cause of 

the corrections is simply the change of the kink shape which is of order 1/C 2 . 



VI. CONCLUSIONS 

We discussed the issue of soliton-assisted energy transport in spatially extended systems 
with external driving and damping. Considering the topological soliton (kink) motion in 
the well-known sin-Gordon model, we showed, that the conventional description of energy 
transport based on the total field momentum does not provide one with the correct value of 
energy flux in the system. In particular, in the case of a directed soliton motion under the 
influence of time-periodic external forces with zero mean, the averaged value of the total field 



momentum is known to be strictly zero 
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21] , while the energy transport associated 



with the soliton motion is obviously non-zero. We identified a new energy pathway - the 
exchange current - which is entirely mediated by the spatial and temporal inhomogeneity 
of the system state. Even for the case of a DC external force, the exchange current is found 
to be small but nonzero. Combining both DC and AC driving we obtain situations when 
the total field momentum is nonzero but the kink does not move on average or moves even 
in the direction opposite to the field momentum. 

The approach to the energy transport in spatially extended systems which is based solely 
on the consideration of the total field momentum reduces the kink motion to that of a point 
particle. However, the soliton motion in presence of AC driving forces is always accompanied 
by the excitation of internal modes [32], which makes the dynamical properties of the moving 
soliton essentially different from those of a point particle. Although we stress here that even 
in cases where no shape modes are excited the total field momentum is not sufficient to 
obtain the total energy current, see Sec. IIII1 It is the spatial extent of the soliton which 
makes the difference for a damped and driven moving soliton. 

The relative contributions of the two current components, the internal and the exchange 
one, can change when considering spatially discrete systems, see Sec. El In particular, for 
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the soliton ratchet case, discreteness induced corrections to the internal current invalidate 
Eq. (0), so that J 1 also contributes to the total energy flow. This explains earlier obtained 
numerical results from Ref. Q]. 

Finally, we would like to note, that our results are also instructive for the general case of 



spatially extendec 
see e.g. Refs. 
Eq 
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systems coupled to external driving forces or other degrees of freedom, 
351 ] . We also mention that the damped and driven sin-Gordon system in 
as well as its discretized version in Eq. ([53)1 . are relevant physical models of annular 
and coupled j^] Josephson junction oscillators, respectively. Therefore, an intriguing 
question rises about the possibility of detection of the new exchange current mechanism on 
the basis of available experimental data (current-voltage characteristics, spectra of emitted 
radiation, etc). Indeed, since directed soliton motion is observed in such experiments, the 
internal and exchange currents have a clear physical meaning, though probably not easy to 



measure. 



APPENDIX 

Using Eq. (|45jl . let us rewrite the expression (|%7j) for the energy current J in the following 
form: 

J = tL f 1 dt [ L dx f dx' [p(x', t) - p{x', t + T)} . (59) 
1 Jo Jo Jo 

According to the multiplicative integration rule |23| 

dx dx' [p(x' , t) — p[x! , t + T)} = / (L — x) [p(x,t) — p(x,t + T)]dx . (60) 
o Jo Jo 

Since the total energy of the system W is a periodic function with period T, the integral 
J L L [p(x, t) - p(x, t + T)]dx = L \W{t) - W{t + T)\ = 0. Thus, expression §§§ is reduced 
to 

J=7^f dt I x[p(x,t)- p(x,t + T)]dx . (61) 
1 * Jo Jo 



Let us introduce the energy-weighted center of kink 

J L x[p(x,t) - p v (t)]dx 



X(t) 



f L [p(x,t)-pv(t)}dx 



(62) 



where p v (t) = p[ip v (t)]. Then we obtain 

T 2 Jo 



J = ±-J T dt f L [p(x,t)-p v (t)]-\x(t + T)-X(t) 



dx . 



(63) 
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Since 



X(t + T) — X(t) = VT, finally we arrive to the following expression for the total 



energy current: 



J= - ! dt f [p(x,t)-p v (t)]dx (64) 
1 Jo Jo 

In a similar way one can derive the compact expression (j51|) for the exchange current. 
We re-write the expression (|48jl in the following form: 

1 r L 



J 1 



T 2 Jo 



dx / dx' 



2T r t 



T 



dt f <f)(x',t')dt' - f dt f <f)(x',t')dt' + T f <p(x',t)dt 
Jt Jo Jo Jo 



Applying twice the multiplicative integration rule, we obtain 

J E = ^ [ L dx{L - x) [ T dt {T<p(x, t + T)-t [<p(x, t + T)- <f)(x, t)}} 
1 Jo Jo 



(65) 



(66) 



By introducing the center of kink X(t) weighted by function <f>(x,t) in analogy to f!62D . 
and using X(t + T) — X(t) = VT, expression (JBTIj) becomes 



J E = [ T dt r (f)(x - VT, t)-(x- Vt)dx - ^ ( l tdt t <p[ip v {t)]dx . (67) 
T Jo Jo T Jo Jo 



Taking into account that (<ft[(p v (t)])T = 0, we finally obtain expression (J5Tj) for the ex- 
change current. 
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